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INTRODUCTION 


The  response  of  germanium  to  shock  waves  should  be  a theoretician's 
delight.  Because  of  its  extreme  importance  to  the  electronics  industry, 
germanium  (Ge)  has  probably  been  studied  more  than  any  other  solid 
material.  Its  basic  (i.e.,  infinitesimal  disturbance)  electrical, 
mechanical,  and  thermal  properties  are  well  understood  in  terms  of  the 
modem  microscopic  concepts  of  solid  state  physics.  Many  mechanical 
(as  well  as  electromechanical)  shock  wave  observations,  however,  have 
not  yet  been  adequately  explained. 

In  this  report  some  of  the  well  known  infinitesimal  disturbance 
properties  of  Ge  are  applied  to  shock  wave  observations.  Calculations 
are  presented  which  strongly  hint  that  electronic  band  structure  effects 
play  a dominant  role  in  at  least  one  (the  approximately  120  kilobars) 
dynamic  pressure-volume  derivative  discontinuity  observation.  If  the 
interpretation  of  the  calculations  is  correct,  a new  view  immediately 
develops  for  the  Hugoniot  elastic  limit  and  the  onset  of  ductility 
(metallic  plasticity)  in  a single  crystal  semiconductor.  That  new  view 
is  correlated  with  the  semiconductor- to-metal  transition  observations 
which  are  characteristic  of  the  diamond  structure. 

Although  silicon  exhibits  shock  wave  properties  very  similar  to 
Ge,  silicon  (Si)  is  not  treated  in  this  report;  more  calculations  need 
to  be  carried  out  before  the  Si  results  can  be  presented. 

The  calculations  presented  contain  detailed  band  structure 
mathematics.  Because  this  type  of  mathematics  is  not  a common  tool  \ 
of  shock  wave  physics,  mini-summaries  are  presented  (clearly  set  off  \ 
from  the  main  text)  which  allow  the  non-mathematically  inclined  reader 
to  appreciate  the  significance  of  the  calculations. 

EXPERIMENTAL  DATA 

Gust  and  Royce  have  presented  a complete  orientational  set  of 
pressure-volume  curves  for  Ge  to  approximately  200  kilobars.  They 
report  lower  stress  "transitions"  occurring  at  approximately  50  kb 
and  approximately  120  kb.  Their  data,  while  containing  approximately 
15  kb  of  scatter,  are  reproduced  in  table  1.  The  data  are  for  a 
uniaxial  strain  shock  propagating  parallel  to  the  indicated  crystal- 
lographic directions. 


1 


s» 

Table  1.  Lower  stress 

shock  transitions  in  germanium3. 

Uniaxial 

Strain 

Direction 

<100> 

<110> 

<lil> 

First  transition 

58  kb 

48  kb 

47  kb 

Second  transition 

117 

119 

125 

a 

From  W.H.  Gust  and  E.B.  Royce,  J.  Appl.  Phys.  £3,  4437  (1972). 

The  common1*  * interpretation  of  the  above  transitions  has  been  to 
identify  the  first  transition  as  the  Hugoniot  elastic  limit  (HFL) , and 
to  identify  the  second  transition  as  the  semiconductor-metal  (polymorphic) 
structural  phase  transformation  as  reported  by  Bundy4 „ Jacquesson  and 
coworkers8,  while  reporting  essentially  the  table  1 data  for  the  <100> 
and  <111>  directions,  also  report  a third  transition  near  130  kb.  Ad- 
ditionally, Jacquesson, et  al5,have  shown  some  hesitancy  to  identify  the 
various  transitions. 

This  report  tentatively  identifies  the  first  transition  as  a mech- 
anical transition  associated  with  brittle  fracture,  and  the  second 
transition  as  ductile  (metallic)  yielding.  The  third  transition  observed 
by  Jacquesson,  et  al* .may  then  turn  out  to  be  the  structural  phase  tran- 
sition observed  (hydrostatically)  by  Bundy4 . 

Before  commencing  with  the  electronic  band  structure  details,  it 
is  worth  noting  that  room  temperature  monocrystalline  Ge  only  undergoes 
brittle  fracture8 *'* 8 . While  twinning8  may  occur  in  room  temperature  Ge, 
that  ductile  twinning  is  always  associated  with  stress  •joncentration  at 
a crack  tip.  The  point  is  that  the  first  yielding  expected  to  be  ob- 
served at  room  temperature  under  shock  conditions  is  brittle  in  character 
- a point  of  great  significance  explained  later  in  this  report. 
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BAND  STRUCTURE  UNDER  UNIAXIAL  STRAIN 


The  room- temperature  atmospheric  pressure  band  structure  of  Ge  is 
shown  in  exaggerated  form  in  figure  1. 


Figure  1.  Room  temperature  - atmospheric  pressure  band  structure  of 
intrinsic  germanium.  From  W.  Paul  and  H.  Brooks9. 
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In  room  temperature-atmospheric  pressure  intrinsic  germanium  the 
four  equivalent  (<111>,  <lll>,  <111>,  and  <1 1 1 >)  valleys  contain 
occupied  electron  states,  while  occupied  hole  states  exist  at  the 
r2j  valence  bands  (the  J = 3/2  states)  near  k = o. 

The  presence  of  shallow  (with  respect  to  Tj)  donor  states  with  the 
0.67  ev  gap  can  (and  usually  does)  increase  the  number  of  electrons  in 
the  fj  state  by  orders  of  magnitude  as  compared  to  the  intrinsic  case. 


The  variation  of  figure  1 with  strain  is  calculated  below.  For 
point  of  reference,  hydrostatic  strain  is  considered  first. 

The  j^1  originally  degenerate  valley  of  the  it^1  type  (e.g.,  Ti  or 
A, ) conduction  band  undergoes  a change  in  energy  minima,  6E.(i),  given 
by  10  11  J 


J 


where  and  5^ 

d u 


II  ♦ 5Ci)  ACi)  A(i) 

d ~ u j j 


(1) 


are  independent  coefficients  which  are  constant 


for  each  band,  e is  the  strain  tensor,  II  is  the  unit  dyad  (unit 
diagonal  3x3  matrix),  and  At1! is  the  unit  vector  characteristic  of 

the  major  axis  of  the  j ^ valley  of  the  i**1  conduction  band.  Repeated 
indices  are  not  to  be  summed  over,  and  6E^  represents  an  energy  shift 

with  respect  to  the  original  energy  level  (or  equivalently  with  respect 
to  a fixed  valence  band) . Equation  1 holds  for  <111>  or  <100>  type 
valleys,  but  not  for  <110>  type  valleys10. 

Denoting  <1 1 1 > type  valleys  by  the  superscript  1,  and  <100>  type 

te 
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valleys  by  the  superscript  2,  allows  one  to  wri 
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and  for  the  six  <100>  type  valleys  (the  relatively  unoccupied  Aj 
band  in  Ge) 


42)  = - 42)  = , As(2)  = - Aj2)  = t2  , A^2)®  - aJ2),  g3  f 


where  denotes  a unit  vector  along  the  i^  crystal  cube  axis  (e.g. 

6*  is  the  unit  vector  parallel  to  the  <010>  cube  edge). 

The  two  J = 3/2  bands  which  are  degenerate  (fig.  1)  at  k = o 
undergo  a strain  induced  energy  change. 


. at  k = o given11  by 

y 

= a l : e + [(Ell  - e22)J  + cp  ] + d2  (e2y  + cp),  ( 

where  a,  b,  and  d are  Material  constants,  and  cp  denotes  cyclic 
permutation. 

Hydrostatic  Compression 

In  keeping  with  the  usual  shock  wave  convention,  this  report 
takes  strain  as  positive  in  compression.  Thus  for  hydrostatic 
compression, 

100  \ 

& = e 010  (2 

°\ 001 / 

Substituting  equation  3 into  equations  1 and  2 gives 

dE^  = 3 5^  + , all  i (; 

i d 3 ° 


6E^  = 3 ^ e , all  i 

i Id  3 ° 


6 E^  = 3 a e . 


Thus  none  of  the  originally  degenerate  bands  split  in  the  hydro- 
static compression  case  (as  expected,  because  hydrostatic  compression 
maintains  the  original  crystal  symmetry) . 


5 values  exist  as  independent  data  (from  acoustoelectric1 2 and 
u , , 

piezoresistance  experiments),  while  S values  do  not  exist  independ- 
ently. Optical  band  gap  transition  experiments9,11  measure 

[dE^1,2^  - 6E^ 

Li  v • 

From  the  Paul  and  Brooks*  data^ 


5.0  ev. 


(7) 


J 2) 

=d 


-1.5  ev. 


(8) 


Error  values  have  been  omitted  from  equations  7 and  8.  The  effect  of 
such  experimental  errors,  although  not  important  to  the  concepts  being 
presented  here,  is  treated  in  the  analysis  of  <I10>  uniaxial  strain 
induced  energy  shifts. 


In  this  report  all  calculations  are  carried  out  for  stress 
(pressure)  values  of  50  kb.  As  only  essentially  linear  theory  (see 
the  < 1 10 > strain  summary)  is  being  considered,  simple  scaling  yields 
predictions  for  any  stress  value.  Employing  a bulk  modulus  value, 

K = 7.50  x 1011  cgs,  and 


P 


gives  for  50  kb 


AV 


= 3 e 

o 


e {50}  = 2.22  x 10'2 
o 


(9) 

(10) 


+In  converting  the  Paul  and  Brooks9  data  to  a strain  basis  a bulk 
modulus  value  of  K = 7.50  x 1011  cgs  has  been  used.  The  Paul  and 
Brooks  data  are  used  here  because  it  is  based  on  room  temperature, 
and  equations  7 and  8 are  known11  to  exhibit  a strong  temperature 
dependence.  The  Balslev11  data  are  for  80  K. 


Substituting  equations  7-10  into  equations  4-6  yields  the  energy  level 
shifts  illustrated  in  figure  2. 
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Figure  2.  Energy  level  shifts  for  50  kb  hydrostatic  compression. 
Energies  are  in  electron  volts. 

The  <11 1>  and  <100>  band  crossing  at  approximately  35  kb  has  been 
verified  experimentally14  (and  has  been  well  understood  for  many 
years) . The  <100>  band  is  predicted  to  collapse  into  the  valence 
band  at  approximately 

I 0.80  \ 

I 0. 10  X 50  = 400  kb. 

Hydrostatic  electrical  resistivity  (room  temperature)  experiments4 
show  a semiconductor-to-metallic  transition  at  approximately  115  kb. 

Hydrostatic  Strain  Summary 

The  experimentally  observed  115  kb  hydrostatic  semiconductor- 
metal  transition  is  in  no  way  correlated  with  any  band  crossings. 

Thus  that  hydrostatic  transition  must  be  associated  with  a polymorphic 
(structural)  phase  transition  wherein  the  new  phase  has  a collapsed 
conduction  band-valence  band  structure.  Furthermore  the  near  room- 
temperature  115  kb  transition  is  quite  sluggish  (ref.  4,  fig.  5), 
casting  some  doubt  on  the  ability  to  observe  it  under  dynamic  shock 
wave  conditions . 


The  hydrostatic  experiments4-  14  do  not  show  any  evidence  of 
the  approximately  50  kb  uniaxial  strain  transition  (table  1)  observed 
under  shock  wave  conditions.  This  is  a direct  indication  that  shear 
(which  is  not  present  hydrostatically)  is  important  for  the  50  kb  shock 
transition. 

<1 1 1 > Uniaxial  Strain 

The  strain  tensor  for  < 1 1 1 > uniaxial  strain  is 

ei  1 1 


e 


(ID 


with  exli  being  the  uniaxial  strain  in  the  <1 1 1 > direction.  Substituting 
equation  11  into  equations  1 and  2 gives 
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The  procedure  now  is  to  evaluate 
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and  then  to  calculate  the  <1 1 1>  strain- induced  hand  splitting  with 
respect  to  the  average  band  motion  of  equation  16. 


If  Tm  is  the  <111>  directed  stress,  then15 


‘in 


an  + 2<Jn  > 


(1«) 


where  o is  the  stress  measured  in  the  cube  axis  coordinate  system, 
ij 

The  a. . are  related  to  the  cube  axis  coordinate  strains,  and  thus 
ij 

Cj, j through  equation  11,  by 
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so  that 
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Employing  the  elastic  constant  values 


c.l 

= 1.29 

X 

101* 

cgs 

Cn 

= 4.83 

X 

10” 

cgs 

C44 

= 6.71 

X 

1011 

cgs 

gives  an  em  value  corresponding  to  50  kb  (for  Tm)  of 
et x i {50}  = 3.64  x 10'*  . 


(19) 


(20) 


(21) 


Equations  7 and  8,  along  with  equation  22,  then  suffice  to  determine 
the  average  band  motion  of  equations  16  and  17. 


The  <1 1 1>  band  splitting  is  determined  by  subtracting  equation  16, 
i a = o,  from  equations  12  and  13.  The  result  is 
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with  a = o,  from  equations 
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The  deformation  potential  is  known18  to  exhibit  a temperature 

dependence  (although  apparently  not  as  strong  as  that  of  the  indirect 
band  gap  values  of  equations  7 and  8 ) . In  this  report 

5(1)  = - 16.7  ev  (24) 

u 

(minus  sign  for  strain  positive  in  compression)  is  used.  That  value 
is  based  upon  piezoresistance  measurements  as  a function  of  temperature 
by  Morin,  Geballe,  and  Herring17,  and  should  be  fairly  representative 
of  room  temperature.  The  16.7  ev  value  is  consistent  with  relaxation 
time  experiments1 8 and  is  essentially  equal  to  the  80  K Balslev11  value 
of  (16.2  _+  0.4)  ev. 

The  average  band  motions,  when  combined  with  equations  21-24,  and 
the  splitting  component  of  equation  15  give  the  50  kb  band  structure 
illustrated  in  figure  3.  In  arriving  at  figure  3,  the  Balslev11  value 
of  d = 3.7  ev  has  been  used  for  the  strain- induced  valence  hand  splitting. 
Figure  3 predicts  that  the  conduction  band  collapses  onto  the  valence 
band  at  a <11 1>  uniaxial  strain  corresponding  stress  of 


<11 l>  Strain  Summary 

Fairly  well  established  deformation  potential  data  predict  a 
valence  band  and  <1 1 1 > conduction  band  crossing  (collapse)  at 
approximately  112  kb  - a value  which  is  in  very  good  agreement  with 
the  observed  125  kb  shock  transition  of  table  1. 

The  hydrostatic  stress  corresponding  to  112  kb  is  simply 
given  by 


P„,  = * Ei 1 1 (112)  , (25) 

or 

P.  12  = 7.5  x 10"  x H2L_  x 3.64  x 10'2  , 

50 


PUJ  =61  kb.  (26) 

Equation  25  neglects  the  effects  of  any  type  of  "yielding  which  may 
be  associated  with  the  observed  50  kb  shock  wave  transition  (table  1). 
Such  yielding  serves  as  a strain  (stress)  release  mechanism  leading 
to  a hydrostatic  stress  less  than  that  of  equation  26.  In  any  event, 

61  kb  is  considerably  less  than  the  experimentally  observed4'  14  hydro- 
static pressure  of  115  kb  necessary  for  the  hydrostatic  semiconductor- 
to-metal  transition.  Furthermore,  because  the  hydrostatically  observed 
transition  is  slow  (i.e.  "sluggish")  it  might  not  be  expected  to  be 
observable  in  microsecond  observation  time  shock  wave  experiments. 

The  numbers  in  figure  3 suggest  the  following  explanation  for 
the  shock  wave  observations: 

1.  Some  type  of  yielding  transition  (which  is  TENTATIVELY 
identified  as  brittle  fracture)  at  50  kb,  followed  by 

2.  Ductile  (metallic)  yielding  which  is  enabled  by  the 
uniaxial  strain  induced  semiconductor- to-metal  (via  the  band  gap 
going  to  zero)  transition  at  approximately  112  kb  of  <1 1 1 > uniaxial 
equivalent  strain. 

In  the  shock  wave  experiment  the  metallic  state  is  achieved  before 
the  polymorphic  transition  occurs.  In  the  hydrostatic  experiment  a 
polymorphic  transition  occurs  prior  to  conduction  band  collapse,  and 
the  band  structure  associated  with  the  new  phase  is  undoubtedly  metallic. 
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The  consequences  of  the  observed  50  kb  transition  have  been 
neglected  in  the  scaling  to  the  120  kb  region.  The  apparent  elastic 
constants  can  undergo  large  changes  when  a mechanical  transition  is 
traversed.  That  neglect  represents  a potentially  important  difficulty. 
Because,  however,  plasticity  does  not  result  in  a net  volume  change, 
plasticity  associated  with  the  50  kb  transition  only  indirectly  affects 
the  strain  to  be  associated  with  the  deformation  potential  uniaxial 
strain  calculations.  This  question  of  elastic-plastic  is  considered 
further  in  the  <100>  uniaxial  strain  summary  of  this  report  (where 
the  calculations  are  the  least  complicated).  While  much  more  work 
needs  to  be  done  on  this  elastic-plastic  question,  it  appears  that 
the  50  kb  transition  can  be  neglected  in  scaling  at  least  to  zeroth 
order . 

<100>  Uniaxial  Strain 

Here  the  strain,  again  referred  to  a cube  axis  coordinate  system, 
is 


lioo\ 

t * Eiool  000  I • 

\ooo  / 

And  upon  again  employing  equations  1 and  2 


i 

j 


=(D 


1 


S(2)  + SC2) 
“d 


u 


f"  )ei##,  311  i* 

loo  » j = 2 • 


:C2) 


E100 


, k = 3,4, 5,6. 


(27) 


(2R) 

(29) 

(30) 


In  equations  29  and  30,  the  j and  k values  refer  to  the  individual 
<100>  type  valley  minima  discussed  previously. 

<5Ev~  * a e100  ± |b  e100  | (31) 

Fot  <100>  uniaxial  strain  only  the  <100>  minima  and  the  valence 
band  are  split.  The  valence  band  splitting  is  given  by  the  second 
term  on  the  right  of  equation  31,  and  the  <100>  splitting  with  respect 
to  the  mean  energy  is  given  by 


Eioo  » j = 1*2. 


(32) 
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I 


I 


4eC2)  < . I =<2>  e,„  , k ■ 3.4, 5,6.  (33) 

k 3 u 

The  value  of  for  Ge  has  not  been  determined  experimentally. 

There  does,  however,  appear  to  be  great  similarity  between  the 
deformation  potential  coefficients  for  Ge  and  Si.  The  hydrostatic 
motion  of  the  <100>  indirect  band  gap  in  Ge  is  equal*  to  that  in  Si 
(where  the  <100>  minima  are  occupied) . Thus 

S(2)  = - 8.6  ev,  (3 4) 

u 

(minus  sign  for  strain  positive  in  compression)  is  taken  from  Si 
experimental11  values  and  used  here.  The  Balslev11  value  of  b * 1.8 
ev  is  used  in  equation  31. 

The  strain  e100  corresponding  to  50  kb  is  given  by 


e100  {50}  = 3.88  x 10- 


The  results  of  completing  the  calculations  indicated  above  are 
illustrated  in  figure  4. 


Figure  4.  Energy  level  shifts  for  <100>  uniaxial  strain  (correspond- 
ing stress  of  50  kb). 
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■ 


I 


4 


j 

\ 
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<100>  Strain  Summary 

Although  the  figure  4 results  are  slightly  suspect  because  of 

(2) 

the  use  of  the  S silicon  value,  these  results  are  entirely  con- 
u 

sistent  with  the  overall  point  of  view  being  presented  in  this  report. 
Two  of  the  <100>  conduction  band  minimas  are  predicted  to  collapse 
onto  the  upward  moving  split  J * 3/2  valence  band  at  a <100>  uniaxial 
stress  level  corresponding  to 


x 50  kb  =*  102  kb. 

I (0.80  - 0.41)  J 

The  above  103  kb  value  is  consistent  with  the  table  1 experi- 
mental value  of  117  kb  for  the  <100>  observed  shock  transition.  It 
is  also  consistent  with  the  observed  ordering  of  the  <100>  and  <11 1> 
shock  transitions  in  that  the  numbers  of  this  report  predict  a con- 
duction band  collapse  for  <100>  uniaxial  strain  at  a lower  stress 
than  the  112  kb  prediction  for  the  <111>  uniaxial  strain  case.  It  is 
necessary,  however,  to  urge  caution  because  of  the  somewhat  guessed 

value  for  5^  . 
u 


As  mentioned  in  the  <1 1 1>  uniaxial  strain  summary,  this  report 
neglects  elastic-plastic  effects  associated  with  the  50  kb  transition. 

Above  the  50  kb  transition  the  total  uniaxial  strain  has  elastic 
and  plastic  components.  If  "e"  denotes  elastic  and  "p"  plastic,  then 
the  total  strain  components  may  be  written  as 


e P 

♦ e^  = o (uniaxial  condition), 

P P 

en  * 2e^  * ° (no  P^stic  volume  change), 

where  11  (parallel)  denotes  the  (<100>,  <100>)  strain  component.  L 
denotes  a perpendicular  (e.g.  (<010>,  <010>))  component. 
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Only  the  elastic  strains  are  important  for  volume  changes 
associated  with  the  deformation  potential  calculations.  Thus 
equation  27  should  be  replaced  by 


above  the  50  kb  transition. 


e®  o o 


can  be  rewritten  as 


/ e 

e 

' eu 

" el 

0 

o 

+ 

o 

o 

o 

L 

o 

/ 

o 

Thus  the  ee  for  deformation  potential  considerations  may  be  viewed 
as  being  composed  of  a hydrostatic  elastic  compression  plus  a uniaxial 
strain  smaller  than  that  predicted  by  equation  36  if  e®  is  positive. 


As  the  hydrostatic  and  uniaxial  parts  of  iee  are  additive,  the 
results  remain  virtually  unchanged  (other  than  for  the  effect  of  the 
e^  components).  Nevertheless,  elastic-plastic  effects  represent  an 

important  area  of  future  work  (especially  in  the  <1 1 1>  and  <110> 
uniaxial  strain  geometries) . 


<110>  Uniaxial  Strain 


The  cube  axis  referred  strain  is  given  by 


1 

Substituting  equation 

37  into  equations  1 and  2 gives 

j 

(SE^  = 
i 

+ 2 Bd)  I - i - l 3 

“d  3 “u  110’  " ’ ' 

(38) 

• 

5e!»  = 

J 

SCU  e 1=12 

-d  E110*  J 

(39) 

s -(2>  1 

6E(2^  = 
i 

^ j e,„.  i = 1. 2,5,4. 

(40) 

SE^  = 

J 

E110»  j = 5,6, 

(41) 

where  the  i and  j values  refer  to  the  specific  valley  minima  identified 
earlier  (following  equation  2). 


4E^±)  = a eixo  1 + 2d2  fill 

V 2 2 


4EW  . 

el  10 

a e,,,  + 

Jb2  + 2d2. 

(4.3) 

V 

2 

V 

The  <111> 

band  splitting  with 

respect  to  the  mean 

is  given  by 

Apf1)  = _1_ 

1 3 

e i - 1 3 

“u  eno*  1 " A»-'» 

(44) 

- - 1 
j 

ell0,  j = 2,4. 

(4  5) 

The  <100>  band  splitting  with  respect  to  the  mean  is 


Al'f2)  = I S(2)  e110,  i = 1,2,3, 4, 
l 6 u 110 


acC2)  _ 1 „( 2)  . . , 

AF . £ Cj  jo*  J “ 5|  6 . 
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The  uniaxial  strain  corresponding  to  the  50  kb  shock  wave  experiment 
is  determined18  by 


Tim  = °n  + °ii 


°n  = (cn  + cm) 


» °12  = °** 


where  o is  the  stress  tensor  measured  in  the  cube  axis  coordinate 

ij 

system.  Combining  equations  48  and  49  yields 


cn  + C12  + 2 c4« 


e110  {50}  = 3.21  x 10"*.  (SI) 

Combining  the  above  <110>  strain  equations  with  the  deformation 
potential  values  expressed  previously  results  in  figure  5. 

<110>  Strain  Summary 

Figure  5 predicts  that 

1.  Two  <1 1 1>  type  valleys  will  collapse  onto  a split-off 
valence  band  at  a <110>  uniaxial  strain  equivalent  stress  of 


0.67'-  0.-56 1 X 5°l  kb  = 30S  kb. 

2.  Four  <100>  type  valleys  will  collapse  onto  a split-off 
valence  band  at  a <110>  uniaxial  strain  equivalent  stress  of 

I °-80  I x 50/  kb  = 211  kb. 

) 0.80-0.19/  ) 

It  is  thus  seen  that  the  <110>  uniaxial  strain  prediction  is  too  high 
(211  kb  vs  the  125  kb  experimental  value  of  table  1).  Assuming  that 
the  model  proposed  in  this  report  is  correct,  the  reason  for  the 
large  <110>  uniaxial  strain  experimental-theoretical  discrepancy  is 
not  clear. 
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As  previously  mentioned,  there  is  some  uncertainty  in  the 
various  deformation  potential  coefficients.  The  theory  presented  here 
is  exactly  linear  in  the  <1 1 1 > and  <100>  uniaxial  strain  cases,  and 
essentially  linear^  for  small  errors  in  the  <1 1 0>  uniaxial  strain  case. 
The  effect  of  coefficient  errors  is  thus  to  raise  (or  lower)  the  con- 
duction band-valence  band  collapse  stress  equally  (e.g.  by  10%  for  a 
10",  uniform  coefficient  error)  for  the  three  uniaxial  strain  directions. 

' 

Only  in  the  circumstance  where  the  errors  in  the  b and  d 
valence  band  coefficients  are  opposite  in  sign  can  a non-uniform 
effect  arise.  In  such  a circumstance,  it  is  possible  to  have  the 
valence  splitting  contribution  to  the  211  kb  value  remain  unchanged 
while  valence  contributions  to  the  <1 1 1>  and  <100>  strain  collapse 
cases  could  be  increased  or  decreased  (depending  upon  the  sign  of  the 
error  in  the  relevant  coefficient).  Since,  however,  the  valence  band- 

j splitting  does  not  represent  the  major  contribution  to  the  energy  gap 

closure  for  any  of  the  three  uniaxial  strain  configurations,  it  is 
difficult  to  see  how  the  relative  ordering  (and  approximate  magnitudes) 
of  our  <1 1 1>,  <100>,  and  <110>  collapse  predictions  can  be  greatly 
affected  by  errors  in  the  deformation  potential  coefficients. 

UNIAXIAL  STRAIN  DUCTILE  YIELDING 

I 

The  "octahedral"  shearing  stress1®  may  be  written  in  the  form 

j - 02  2 ) 2 + CP>  + b {°l2+  CP>  * (52) 

| t is  the  "resolved"  shear  stress  in  the  < 1 1 1 > plane, 

oct 

Alexander  and  Haasen7  quote  25  kb  as  the  ( < 1 1 1 > plane)  shear  stress 
necessary  to  create  a dislocation  in  Ge.  The  procedure  is  thus  to  set 

t = 25  kb  and  to  calculate  the  equivalent  uniaxial  strain  correspond- 
oct 

ing  stress. 


^Substitute  b'  = b + e and  d'  = d + e into  equation  43,  and  expand 

b d 

the  square  root.  The  errors  are  then  seen  to  enter  the  final  result 
linearly,  and  to  have  a maximum  effect  on  the  < 1 10 > uniaxial  strain 


20 


<100>  Uniaxial  Strain 


Here,  using  equation  27,  it  is  easily  seen  that 

CTn  cn  eioo  » °2  2 = ajj  = C12  eioo  » C***-^) 


so  that 


r . 'll  (C  - c ) e 

oct  3 ^11  loo  * 


or 


JT 


oct 


C1  2 

1 - r i °n 

'"X  1 


(541 


(55) 


Setting  r = 25  kb  and,  using  the  previously  listed  elastic  constant 
oct 

values,  gives 


hi  = • 83  kb 

0.30 


(56) 


<lll>  Uniaxial  Strain 

From  equations  18  and  19,  and 

e 


t 


‘i  i l 


= c 


xy 


44 


(57) 


it  is  easily  seen  that 


3/2 


Toct  ’ 3 C,‘  Cl" 


(58) 


+The  factor  of  2 appears  in  equation  57  because  Chou  and  Pagano15  use 
a shear  strain  twice  that  associated  with  equations  11  and  37. 
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ant’  finally  employing  equation  20  gives 


' I 


3/2 

_ 2 C44 


oct  Cjj  + 2 c12  + 4 c44 


Tm  > 


(59) 


where  Tj j , is  the  <111>  directed  stress  corresponding  to  <1 1 1 > uniaxial 


stress . 


Thus  upon  using  the  previously  listed  elastic  constant  values  and 

(60) 


t =25  kb.  one  sees  that 
oct 


S/2/o.67ll 

Vt = 2 (or | 


or 


tj  x j = 2.60  t = 65  kb. 


oct 


<lin>  Uniaxial  Strain 

From  equations  37,  48,  49,  and  52 


aU  - °22  = tCll  + CI2) 


ei  10 


"33  = 2 C12 


110 


°12  = C44  E110 


'23 


°S1  = ° » 


°1  1 0 


(Cli  - Cl2)' 


oct  3 V 

or  upon  using  equation  50 

T 


+ 6 c 


4 4 > 


'110 


N 


(Cll  * C12  + 2 C44)toct 


(Cn  " C12)1 


♦ 6 cJ 

44 


(61) 


(62) 


(63) 


(64) 


22 


1X0 


3.113 

1.74 


oct 


["ijo  = 67  kb. 


where  x = 25  kb  has  again  been  used, 
oct 


(65) 

(**) 


There  are  two  points  to  be  made  concerning  the  numerical  results 
displayed  in  equations  56,  61,  and  66: 

1.  The  <111>  plane  shear  stress  necessary  for  dislocation 
generation  corresponds  to  a uniaxial  strain  equivalent  stress  at  least 
30%  higher  than  the  observed  50  kb  first  shock  wave  transition  (table  1) . 

2.  Equations  56,  61,  and  66  are  consistent  with  conduction  band- 
valence  band  collapse  being  the  switch  which  allows  the  lower  stress 
(65,  67,  and  83  kb)  dislocation  generation  to  proceed  in  a ductile 
medium. 


CONCLUSIONS 

A rough  equivalency  has  been  demonstrated  between  the  uniaxial 
strain  corresponding  stress  for  conduction  band-valence  band  collapse 
and  the  approximately  120  kb  observed  shock  wave  transition.  If 
difficulty  with  the  calculations  exists,  it  is  to  be  found  in  the 
high  (211  kb)  claculated  stress  associated  with  the  <110>  uniaxial 
strain  case. 

The  TENTATIVE  conclusion,  based  upon  the  observation  that  brittle 
fracture  precedes  ductile  yielding  in  room  temperature  Ge,  is  that  the 
50  kb  shock  wave  transition  is  brittle  in  character.  That  brittle 
fracture  then  is  followed  by  a conduction  band- valence  band  collapse 
at  approximately  120  kb  (211  kb  for  <110>  uniaxial  strain)  which  acts 
as  an  enabling  switch  for  ductile  (metallic)  yielding. 

As  mentioned  previously,  the  presence  of  brittle  fracture  at  50  kb 
must  certainly  (in  that  it  acts  as  a stress  relief  mechanism)  alteT  the 
strain  matrices  from  those  used  in  uniaxial  strain  calculations  pre- 
sented here.  While  it  has  been  shown  that  the  50  kb  elastic-plastic 
effects  cancel  to  zeroth  order,  those  effects  represent  an  area  where 
more  theoretical  work  must  be  performed.  It  is,  however,  worth  noting 
that  perhaps  a fundamental  and  important  difference  exists  between  the 
consequences  of  brittle  and  ductile  yielding;  brittle  fracture  can  be 
thought  of  as  resulting  in  pieces  (cleaved)  of  the  original  material. 


I 
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with  each  piece  being  essentially  only  elastically  deformed.  Ductile 
yielding  should  be  a more  extensive  phenomenon  allowing  for  the  concept 
of  a uniform,  plastic-strain  contribution. 

To  the  extent  that  a dislocation  phenomenon  is  involved  in  ductile 
yielding  and  brittle  fracture  (e.g.  via  dislocation  twinning8)  it  is 
expected  that  the  impurity  content  will  affect  the  magnitude  of  the 
yielding  stresses;  dislocation  velocity  is  known7  to  depend  upon  both 
the  number  density  and  the  type  of  impurity.  Because  the  character  of 
an  impurity  can  easily  be  changed  by  the  strain  induced  shifting  of  an 
energy  band  (e.g.,  by  going  from  a neutral  to  an  ionized  donor  state  as 
a conduction  band  minima  decreases  in  energy)  it  is  expected  that  impurity 
content  will  also  play  a role  in  the  stress  magnitude  associated  with 
the  ductile  and  brittle  yielding. 

Perhaps  the  impurity  phenomenon  could  selectively  reduce  the  211  kb 
value  for  the  conduction  band-valence  band  collapse  under  <110>  uni- 
axial strain.  Unfortunately  it  is  not  now  known  how  to  include  such 
impurity  effects  in  the  type  of  calculations  presented  in  this  report. 
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